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Abstract. According to Landauer’s principle, a minimum amount of energy 
proportional to temperature must be dissipated during the erasure of a classical 
bit of information compensating the entropy loss, thereby linking the information 
and thermodynamics. Here we show that the Landauer limit of energy dissipation 
is achievable in a shape-anisotropic single-domain magnetostrictive nanomagnet 
having two mutually anti-parallel degenerate magnetization states that store a bit of 
information. We model the magnetization dynamics using stochastic Landau-Lifshitz- 
Gilbert equation in the presence of thermal fluctuations and show that on average the 
Landauer bound is satisfied, i.e., it accords to the generalized Landauer’s principle for 
small systems with stochastic fluctuations. 
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l. Introduction 

Computers are physical systems and a computation or processing of information whether 
performed in an electronic machinery or in a biological system is subjected to physical 
principle, e.g., the second law of thermodynamics. According to Landauer’s principle [1], 
in a classical system with two degenerate ground states, it necessarily dissipates an 
energy of amount kTln{2) (~ 3 x 10“^^ joule at T=300 K, k is the Boltzmann constant, 
and T is temperature) to erase a bit of information (logically irreversible computation 
contrary to reversible computation [21 a a El El d), which compensates for the entropy 
loss. In fact, the intimate relationship between information and entropy dates back over 
a century ago when Maxwell introduced his famous and controversial demon iiigiiini 

m, which was later exorcized by Bennett [121 ESI EH EH ESI EH Ell ESI [20]. However, 
the possibility of stochastic violation of Landauer’s principle cannot be discarded in 
small systems that are prone to thermal fluctuations [211 1221123l 1211125l [2S|- Indeed 
the stochastic violations of second law have been experimentally observed mm, 
however, on average the second law is still safeguarded [29]. There are controversies 
whether Landauer’s bound holds for asymmetric memories [HO] Ell 132] and in quantum 
regime [331131], although the second law of thermodynamics remains intact. 

Recently, Landauer limit of energy dissipation is experimentally demonstrated for 
a colloidal particle with linear motion trapped in a double-well potential landscape 
utilizing a laser beam [35] • The heat dissipated during the erasure process (i.e., while 
resetting the particle in one of the two wells) is measured and shown to be approaching 
the Landauer’s limit. However, the Landauer’s limit is not experimented for a rotational 
body like the magnetization in a nanomagnet [36l EH EH [39l EOl IH]. We consider the 
complete three-dimensional potential landscape and full three-dimensional motion of the 
magnetization, unlike assuming an overdamped particle with linear motion [231 [25] and 
show that the magnetization may deflect out of magnet’s plane during its dynamical 
motion. Consideration of even a small out-of-plane excursion of magnetization plays a 
crucial role in shaping the magnetization dynamics while achieving the Landauer limit 
of energy dissipation. Here, we particularly consider a single-domain magnetostrictive 
nanomagnet and show that Landauer limit of energy dissipation is achievable in such 
a magnetostrictive particle. Considering that the excessive energy dissipation during 
switching of bits has been the bottleneck behind utilizing traditional charge-based 
transistor electronics further [3211331133], magnetostrictive nanomagnets in multiferroic 
heterostructures have potential to become the staple of modern non-volatile memory 
and logic systems [35ll36ll37ll3Ell39lEQl[5llE2lE3lE3lE5lE6lE7lEH]. We have 
solved stochastic Landau-Lifshitz-Gilbert (LLG) equation [59l EOl El] in the presence of 
thermal fluctuations to theoretically analyze the aspects of magnetization dynamics and 
determine the energy dissipation during the erasure of a bit of information. Simulation 
results also show that stochastic violation of the Landauer bound is possible, nonetheless 
the bound is respected on average. 
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Figure 1. (a) A shape-anisotropic single-domain magnetostrictive nanomagnet and 
axis assignment, (b) The time-cycle of uniaxial stress (along z-direction) and magnetic 
field, H (along 0 = 0° direction) applied on the nanomagnet, (c) The in-plane {y-z 
plane, (p = ±90°) potential landscapes of the nanomagnet at different instants of 
time, (d) The complete three-dimensional potential landscapes of the nanomagnet at 
different instants of time. The solid lines on the landscapes correspond to the in-plane 
potential profiles, which are shown in (c). Consideration of these out-of-plane potential 
landscapes plays an important role in shaping the magnetization dynamics discussed 
later on. At time ti, magnetization can reside in either of the wells (0 = 0° or 0 = 180°) 
with 50% probability each, however, at time t 2 , it reaches at 0 = 90°, i.e., the bit is 
erased. At time and onwards, magnetization reaches at 0 ~ 0° deterministically. 


2. Model 

We consider that the single-domain magnetostrictive nanomagnet is shaped like an 
elliptical cylinder with its elliptical cross-section lying on the y-z plane; the major axis 
and minor axis are aligned along the 2 ;- and ^/-direction, respectively (see Fig. [1^). In 
standard spherical coordinate system, 6 is the polar angle and 0 is the azimnthal angle. 
Any deviation of magnetization from magnet’s plane {y-z plane, 0 = ±90°) is termed 
as ont-of-plane excnrsion. The dimensions of the major axis, the minor axis, and the 
thickness are a, b, and I, respectively. So the magnet’s volnme is = {7r/4:)abl. Dne 
to this shape anisotropy, the two degenerate states 6* = 0° and 9 = 180° along the 
z-direction (easy axis) can store a bit of information. The y-axis is the in-plane hard 
axis and the a:-axis is the ont-of-plane hard axis. Since I <C b, the ont-of-plane hard axis 
is mnch harder than the in-plane hard axis. 
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Figure [T]d shows the time-cycle of uniaxial stress (along z-direction) and a 
magnetic held H (along 6 = 0° direction) applied on the magnetostrictive nanomagnet. 
The adiabatic pulses should be slow enough that magnetization follows its potential 
landscape and the quasistatic assumption is valid. We can write the total energy of the 
nanomagnet per unit volume as the sum of three energies - the shape anisotropy energy, 
anisotropies induced due to applied stress and asymmetry-making magnetic held (note 
that it is possible to replace the magnetic held by intrinsic interface coupling between 
polarization and magnetization in multiferroic heterostructures [5^154] ) - as follows 133: 

1 4^1 .®s/iape(0) sifl 6 B COS 6 COS 6, (1) 

where 

Bghape{(l)) = (1/2) M[Hk +HdCOS^cj)], (2) 

Bgtress{t) = {3/2) Xga{t), (3) 

Basymmit) = MH{t), (4) 

M = fioMg, po is permeability of free space, Mg is the saturation magnetization, 
Hk = {Ndy — Ndz) Mg is the Stoner-Wohlfarth switching held [36], Hd = {Ndx — Ndy) Mg 
is the demagnetization held [37], Ndm is the m*^ (m = x,y,z) component of the 

demagnetization factor [40] {Ndx S> Ndy > Ndz), (3/2)A., is the magnetostriction 

coefficient of the magnetostrictive nanomagnet 133. a = Te is the stress, Y is the 
Young’s modulus, e is the strain that generates the stress, and H{t) is the magnetic 
held at time t. 

Note that a negative Aga product will favor alignment of the magnetization along 
the minor axis {y-axis) and according to our convention, in a material having positive As, 
a sufficiently high compressive stress that can overcome the in-plane shape-anisotropic 
barrier will rotate the magnetization toward the in-plane hard axis {9 = 90°, 0 = ±90°). 
This is depicted in Figs. and Hji that magnetization initially starting from any of the 
two degenerate stable states {6 = 0°, state 0 and 6 = 180°, state 1 with probabilities, 
Po = Pi = 0.5) comes to the in-plane hard axis at time t 2 - The initial entropy of the 
system at time ti is —kY/nVnlnpn = kln{2). However, when stress is applied between 
time ti and t 2 , the barrier separating the two stable states gets removed and the potential 
landscape becomes monostable in ^-space (at 6 = 90°) with entropy zero. A reduction 
of entropy k ln{2) must be dissipated as heat during this process. The reason that the 
barrier is removed before applying any asymmetry-making held H is to make the erasure 
process independent of barrier height, which determines the hold failure probability and 
also to resist thermal huctuations by making the monostable well deep enough. 

Note that upon reaching 9 = 90°, magnetization may situate in (/-space either at 
(j) = ±90° or at (/ = —90°, which is a consequence of having the planar shape of the 
nanomagnet, however, it should not be related to 6*-space in which the magnetization is 
switched. The (/-space does exist even at 6* ~ 0° or 6* ~ 180° [e.g., see Supplementary 
Fig. SI, where there exist two peaks at (/ = ±90° around 6 ~ 180°]. 
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Between times ^2 and t^, an asymmetry-making field H is applied to 
deterministically rotate the magnetization to state 0 = 0°. The degree of this tilt 
should be sufficient enough to dissuade thermal fluctuations. Afterwards, both the 
stress and H are removed to complete the switching process. Note that the entropy 
of the system at time t 2 and onwards are zero. Hence there is no minimum energy 
dissipation bound when magnetization traverses from time t 2 to as the dissipation 
can be made arbitrarily small. 

The torque acting on the magnetization is derived from the gradient of potential 
landscape. Additionally, there is a random thermal field [h(t) = hi(t) bi, i G (x, y, z)] 
to incorporate thermal fluctuations [61]. Solving the stochastic Landau-Lifshitz- 
Gilbert equation [59l [601 El] of magnetization dynamics analytically (see Supplementary 
material for details), we get the following coupled equations for 6 and <p: 

(l + ^ = ^-^[-(^B{(j),t)sin{2e) - aBasymm{t) sine 

+ Bshape, 4 >{(t))sine {aPe{e, 0, t) + P^{e, 0, f))], (5) 

(l + ^ = ^[25(0,^)COS0 Basymm{t) 

+ aBshapeA^) - [sine]~^ {Pe{e, 0 , t) - aP^{e, 0 , t))] 

{sine 7^ 0), (6) 

where 

B{(j),t) = BshapeiA + Bstressit), A) 

Bs/iape, 0 ( 0 ) = ( 1 / 2 ) Mifrf sm( 20 ), ( 8 ) 

Pg{e,(j),t) = M[hs:{t) cose cos(j) + hy{t) cosesincj) — hz{t) sine], (9) 

P^{e,(j),t) = M[hy{t) coscj) — hx{i) sincf)], (10) 

^ y (* = x,?/,z), (11) 

a is the dimensionless phenomenological Gilbert damping parameter [621 ESI EH ES], 7 
is the gyromagnetic ratio for electrons, 1/Af is proportional to the attempt frequency 
of the thermal held, Atgim = 10^ At is the simulation time-step used, and the quantity 
G(o,i)(t) is a Gaussian distribution with zero mean and unit variance. 

The energy dissipated in the nanomagnet due to Gilbert damping can be expressed 
as Ed = A Pd{t)dt, where r is the time taken during the erasure cycle [i.e., — ti in 

Fig. [11(b)], and Pd{t) is the power dissipated at time t per unit volume given by 

Pi{t) = ite («((). m. t)f, ( 12 ) 

where Te is the torque (see Supplementary material for details) acting on the 
magnetization due to the gradient of potential energy Etotai [Equation ([T])]. Thermal 
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field with mean zero does not cause any net energy dissipation but it causes variability 
in the energy dissipation by scuttling the trajectory of magnetization. 

Assuming magnetization resides in-plane (0 = ±90°) between times ti and ^ 2 , 
(3/2)AsCrmm = (Train IS the stress at which the energy barrier is just 

removed, and Tmax = o'(^ 2 ) = 2amin, we can deduce from Equation ([5]) without 
considering thermal fluctuations (i.e., h = 0 ) that 



(13) 


Note that at 6 *i = 180° and 02 = 90°, the above expression is undehned or it 
would take exceedingly large amount of time to start exactly from 9 = 180° and to 
reach exactly a.t 6 = 90° since the torque acting on the magnetization at these points 
due the gradient of potential landscape vanishes. Fortunately, thermal fluctuations 
can prevent this lockjam in such circumstances. Assuming 6i = 175° and 62 = 95°, 
In(tan92/tan9i) = 4.8725. This closed-form expression gives us an estimate over the 
time required for magnetization to reach 6 ~ 90° depending on the material parameters 
involved. 

It needs mention here that particularly in the presence of thermal fluctuations, 
magnetization may temporarily traverse on higher potential (not only in-plane of the 
nanomagnet but also out-of-plane, i.e., when (j) 7 ^ ±90°) and dissipate energy when it 
comes back to lower potential. A less than one degree of deflection in the out-of-plane 
direction can have immense consequence on magnetization dynamics. The key reason 
behind is that the out-of-plane demagnetization held is about a couple of orders 
of magnitude higher than Hk. When the magnetization dehects out-of-plane due to 
torque exerted on it, an additional torque of comparatively very high magnitudes comes 
into play [see Equation (|H])], which makes the dynamics fast. It is true that the out-of¬ 
plane excursion causes power dissipation but switching also becomes fast, so that the 
net energy dissipation may be smaller compared to the case when it is assumed that 
magnetization is conhned to the magnet’s plane. Since the Landauer limit of energy 
dissipation is very small, we should particularly take into consideration this signihcant 
effect due to out-of-plane excursion of magnetization. 

We will now analyze the magnetization dynamics between times ^2 and t^. 
Assuming (3/2)Xsamax = —MHk and (j) = ±90°, the total energy becomes Etotai = 
— {l/2)MHksin^9 — MH{t)cos6. We see that as H goes from 0 to the minimum 
value of 9 goes from 9 = 90° to 0 = 0° continuously following the expression 


9min{t) = cos ^{H{t)/Hk). 


Due to the smooth transitions of magnetization, with adiabatic pulses between 


times ti and fa, there is no lower bound of energy dissipation in the absence of thermal 


fluctuations. However, the Landauer principle remains intact since this is T = 0 K 
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case and thus the energy dissipation proportional to temperature is also zero. If we 
incorporate random thermal fluctuations at hnite temperature, magnetization will get 
deflected uphill in the potential landscape even with very slow ramp of pulses and will 
incur energy dissipation, which is subjected to Landauer bound. 


3. Results 

We consider the magnetostrictive nanomagnet to be made of polycrystalline Galfenol 
(FeGa), which has the following material properties - magnetostrictive coefficient 
((3/2)A5): +150 x 10“®, saturation magnetization (Mg): 8x10® A/m, Gilbert damping 
parameter (a): 0.025, and Young’s modulus (Y): 140 GPa [66l|67]. The dimensions of 
the nanomagnet is chosen as a = 100 nm, b = 90 nm, and I = 6 nm, which ensures that 
the nanomagnet has a single ferromagnetic domain j39]. With the chosen dimensions, 
the Stoner-Wohlfarth switching held Hk becomes ~0.01 Mg. The values of stress {<Jmax), 
strain (e), and asymmetric held {H) are 60.6 MPa, 433x10“®, and 0.01 T, respectively. 


(a) 


Dynamics of e 


(b) 



Dynamics of ( 



Figure 2. Magnetization dynamics when tn+i — tn = 400 ns (n=l-4). Thermal 
fluctuations are not incorporated during this simulation, however, since exactly at 
9 = 180°, the torque acting on the magnetization is zero, an initial 5 degree of deflection 
for 9 is chosen. The initial value of (p is 90°, choice of the other in-plane angle 4> = 270° 
(—90°) is analogous. Both the cases when magnetization is assumed to be confined on 
magnet’s plane (2D) and that when considering the full three-dimensional potential 
landscape (3D) are plotted, (a) 0-dynamics, and (b) ^-dynamics. 


Figure [2] shows the magnetization dynamics when both the ramp-up and ramp- 
down times of stress and H held are 400 ns (Equation (IT^ gives a value of ~400 ns). 
Both the results - when magnetization is assumed to be restricted on the magnet’s 
plane (2D assumption) and while considering full there-dimensional potential landscape 
- are plotted. Note that even a less than one degree dehection of magnetization out of 
magnet’s plane (see azimuthal angle (/-plots in Fig. |2 ]d) has caused a signihcant change 
in the ^-dynamics due to inherent O-cp coupling in LEG equation. The actual dynamics 
is quite faster than when 2D assumption is enforced. In small systems, the out-of-plane 
demagnetization factor N^x is not that high compared to the y- and ^-component (ratio 
is about 10 times) and since there are torques acting on the magnetization in the cp- 
























direction, magnetization can easily deflect out-of-plane particularly in the presence of 
thermal fluctuations. 


Distribution of energy dissipation, T=300 K 
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Figure 3. Distribution of energy dissipation at room-temperature (300 K) during 
the time interval t 2 — when stress is brought up from zero to a maximum value 
decreasing the entropy of the system by kln{2) = 0.6932fc. A concomitant amount 
of energy must be dissipated according to Landauer’s principle. Note that for a few 
cases, the Landauer bound is violated but the mean of the energy dissipation (0.74 kT) 
exceeds the Landauer bound safeguarding the Landauer’s principle and the second law 
of thermodynamics. A moderately large (1000) number of simulations in the presence 
thermal fluctuations have been performed to generate this distribution. 


Figure [3] shows the distribution of energy dissipation at 300 K during the time 
interval t2 — ti = 100 ns to depict that the generalized Landauer principle of energy 
dissipation is maintained. Stochastic violation of Landauer’s bound 0.6932 kT due 
to thermal fluctuations is observed but the mean energy dissipation does respect the 
Landauer’s bound. In the presence of thermal fluctuations, we note that a slower ramp 
rate of pulse is sufficient than that for 0 K, since thermal fluctuations does a net favor 
to direct magnetization going downhill. The mean energy dissipation during the time 
interval — t2 (when magnetization traverses from 6 = 90° towards 6 = 0°) is 0.07 kT, 
which does not have any bound since there is no entropy loss in the system. 


4. Conclusions 

We have shown that Landauer limit of energy dissipation is achievable in a 
magnetostrictive nanomagnet. With realistic parameters and an well-accepted model, 
i.e., Landau-Lifshitz-Gilbert equation of magnetization dynamics, we have performed 
simulations and the results show that unlike the case of a particle with linear motion, the 
out-of-plane excursion of magnetization plays a crucial role in shaping the magnetization 
dynamics while achieving the Landauer limit. We hope that our finding would stimulate 
experimental efforts to demonstrate that the magnetostrictive nanomagnets are suitable 
for exploring the thermodynamic limit of energy dissipation. On applied perspective, 
such miniscule dissipation can be the basis of ultra-low-energy computing for our future 
information processing systems. Also, this can open up unprecedented applications 
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that need to work with the energy harvested from the environment e.g., monitoring an 
epileptic patient’s brain to report an impending seizure by drawing energy solely from 
the patient’s body. 
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SI. Solution of the stochastic Landau-Lifshitz-Gilbert (LLG) equation 

We consider that the single-domain magnetostrictive nanomagnet is shaped like an 
elliptical cylinder with its elliptical cross-section lying on the y-z plane; the major axis 
and minor axis are aligned along the z- and ^/-direction, respectively (see Fig. la in the 
main Letter). In standard spherical coordinate system, 6 is the polar angle and 0 is 
the azimuthal angle. Any deviation of magnetization from magnet’s plane {y-z plane, 
(j) = ±90°) is termed as out-of-plane excursion. The dimensions of the major axis, the 
minor axis, and the thickness are a, b, and /, respectively. So the magnet’s volume is 
= (7r/4)a6/. Due to this shape anisotropy, the two degenerate states 6 = 0° and 
6 = 180° along the z-direction (easy axis) can store a bit of information. The y-axis 
is the in-plane hard axis and the x-axis is the out-of-plane hard axis. Since I b, the 
out-of-plane hard axis is much harder than the in-plane hard axis. 

We can write the total energy of the nanomagnet per unit volume as the sum of 
three energies - the shape anisotropy energy, anisotropies induced due to applied stress 
and magnetic held (note that we assume a polycrystalline nanomagnet, so that we ignore 
the magnetocrystalline anisotropy energy) - as follows: 


Etotal{6, 0 , t) — Eshape{6, (p) E stress {6, t) ± Easymm{6, t) 

Eshape{4^^ 6 Rstress(^) COS 6 Rasymm(^) COS 6^ (^^) 

where 

Bshape{(p) = {1/2) M[Hk + HdCOS^Cp], (S2) 

Bstress{t) = (3/2)A,cr(t), (S3) 

Basymm{t) = MH{t), (S4) 

M = fioMg, fio is permeability of free sapce. Mg is the saturation magnetization, 
Hk = {Ndy — Ndz) Mg is the Stoner-Wohlfarth switching held [ 91 ], Hd = {Ndx — Ndy) Mg 
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is the demagnetization field Ndra is the (m = x, ?/, z) component of the 

demagnetization factor [32] ^ N^y > N^z), (3/2)As is the magnetostriction 

coefficient of the magnetostrictive nanomagnet m, a{t) is the stress, and H{t) is the 
magnetic field, at time t. 

Note that a negative AsCr(f) product will favor alignment of the magnetization along 
the minor axis (j/-axis) and according to our convention, in a material having positive A^, 
a sufficiently high compressive stress that can overcome the in-plane shape-anisotropic 
barrier will rotate the magnetization toward the in-plane hard axis {9 = 90°, 0 = ih90°). 

The magnetization M of the nanomagnet has a constant magnitude but a variable 
direction, so that we can represent it by a vector of unit norm rim = M/|M| = br where 
br is the unit vector in the radial direction in spherical coordinate system represented 
by (r,6*,0). The other two unit vectors in the spherical coordinate system are denoted 
by bg and b^ for 6 and 0 rotations, respectively. 

The effective field and torque acting on the magnetization per unit volume 
due to gradient of potential landscape can be expressed as Heff = —VEtotai = 
-{dEtotai/dO) be - {l/sin9) {dEtotai/d4>) ^<t> and Te = rim x Heff, respectively. 


Te( 6 ', 0, t) = —[5(0, t) sin{29) + B 

asymm (t) sinO] b^ - Bshape, sinO be, 

(S5) 


where 


Eshapei,4^') T Bgtressii^) 
BshapeA<P) = ( 1 / 2 ) M5rfsm(20). 


(5 6 ) 

(57) 


The effect of random thermal fluctuations is incorporated via a random magnetic 
field h(f) = hxit)e:^^ + hy{t)ey + hz{t)e2, where hi{t) {i = x, y, z) are the three components 
of the random thermal field in Cartesian coordinates. We assume the properties of the 
random field h(t) as described in Ref. [33] . The random thermal field can be written 
as [33] 



(S 8 ) 


where a is the dimensionless phenomenological Gilbert damping parameter, 7 is the 
gyromagnetic ratio for electrons, 1/At is proportional to the attempt frequency of the 
thermal field, Atgim = 10^ At is the simulation time-step used, and the quantity G(o,i)(t) 
is a Gaussian distribution with zero mean and unit variance. 

The thermal field and corresponding torque acting on the magnetization per unit 
volume can be written as Hth( 6 ', 0 , t) = Pe{6, 0 , t) bq -|- P^{6, 0 , t) b^ and Tth( 6 ', 0 , t) = 
Hm X Hth( 6 ', 0 , t), respectively, where 


Fe( 6 *, 0, t) = M[hx{t) cosO coscf) -|- hy{t) cosOsincl) — h^it) sinO], 
P^{9,(f),t) = M[hy{t) cos(j) — hx{t) sm0]. 


(S9) 

(SIO) 
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The magnetization dynamics under the action of the torques TE(t) and Tth(^) is 
described by the stochastic Landau-Lifshitz-Gilbert (LLG) equation [Sl3]-S|5] as follows. 

^ “ i"” liTj “ "m 

Solving the above equation analytically, we get the following coupled equations of 
magnetization dynamics for 6 and 0: 

(l + ^ = ^-^[-(^B{(j),t)sin{2e) - aBasymm{t) sin9 

+ Bshape,4,{^)sin9 + {aPe{9, cj), t) + P^{9, 0, t))], (S12) 

(l + ^ = ^-^[‘^B{(j),t)cOs9 + Basyram{t) 

+ aBshape,4,{4>) - [sin9]~^ {Pe{9, 0, t) - aP^{9, 0, t))] 

{sin9 7 ^ 0), (S13) 

We need to solve the above two coupled equations numerically to track the trajectory 
of magnetization over time, in the presence of thermal fluctuations. 

The energy dissipated in the nanomagnet due to Gilbert damping can be expressed 
as Ed = Jq Pd(t)dt, where r is the switching delay and Pd(t) is the power dissipated at 
time t per unit volume given by 

ft(«) = |Te(»(«). m. (S14) 

Thermal held with mean zero does not cause any net energy dissipation but it causes 
variability in the energy dissipation by scuttling the trajectory of magnetization. 


S2. Distribution of initial orientation of magnetization 


When the magnetization direction is exactly along the easy axis, i.e., sin 0 = 0 (6^ = 0° or 
9 = 180°), the torque acting on the magnetization given by Equation fISSp becomes zero. 
That is why only thermal huctuations can dehect the magnetization vector exactly from 
the easy axis. Gonsider the situation when 9 = 180°. From Equations flS12D and fISlOp . 
we get 


0(t) = tan ^ 


/ ahy{t) + hx{t) \ 
\hy{t) - ah^{t)) ’ 


(S15) 


d9{t) _ _ -\l\{hl{t) + hlit)) _ 

dt ^/{hy{t) - ahs,{t))‘^ + (ahyit) + 

From the Equation fISlOl) . we notice clearly that the thermal huctuations can dehect 
the magnetization exactly from the easy axis since the time rate of change of 9{t) is 
non-zero in the presence of thermal huctuations. Note that d9{t)/dt does not depend 
on the component of the random thermal held along the z-axis, i.e. hz{t), which is 
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a consequence of having z-axis as the easy axis of the nanomagnet. However, once 
the magnetization direction is even slightly deflected from the easy axis, all three 
components of the random thermal held along the x-, y-, and z-direction would come 
into play. 


Distribution of e, T=300K 



(a) 


Distribution of 4., T=300K 



(b) 

Figure SI. Distribution of polar angle 6 and azimuthal angle (p due to room- 
temperature (300 K) thermal fluctuations (a) Distributions of the polar angle 

9. This is a Boltzmann distribution with mean value ~175°, while the most likely 
value is 180°. (b) Distribution of the azimuthal angle p. These are two Gaussian 
distributions with peaks centered at 90° and 270° (or —90°), which depicts that the 
most likely location of the magnetization vector is in the plane of the nanomagnet. 


When neither stress nor asymmetry-making H held is applied on the 
magnetostrictive nanomagnet, magnetization just huctuates around an easy axis 
provided that the shape anisotropy energy barrier is enough high to prevent spontaneous 
reversal of magnetization from one stable state to another in a short period of time. 
We have solved the Equations flS12l) and flSldll while setting a = H = 0 to track the 
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dynamics of magnetization due to thermal fluctuations. So this will yield the distribution 
of the magnetization vector’s initial orientation when stress or H is turned on. Since 
the most probably value of magnetization is along an easy axis, the ^-distribution is 
Boltzmann peaked at 6^ = 0° or 180°. The 0-distribution is Gaussian peaked at 0 = ±90° 
since these positions are minimum energy positions. Taking the initial distribution into 
account does introduce variability in the metrics calculated, however, it does not change 
the mean value signihcantly since the odds are exponentially suppressed. 
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